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Abstract 

A method for the estimation of cross-bispectra of discrete real time 
series is developed. The asymptotic variance properties of the bispectrum 
are reviewed, and a method for the direct estimation of bispectral variance 
is given. The symmetry properties are described which minimize the compu- 
tations necessary to obtain a complete estimate of the cross -bispectrum in 
the right -half-plane. A procedure is given for computing the cross- 
bispectrum by subdividing the domain into rectangular averaging regions 
which help reduce the variance of the estimates and allow easy application of 
the symmetry relationships to minimize the computational effort. As an 
example of the procedure, the cross-bispectrum of a numerically generated, 
exponentially distributed time series is computed and compared with theory. 
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1, Introduction 


Bispectral analysis is currently experiencing a degree of popxdarity with 
applications in a number of diverse fields. The bispectrum was originally pro- 
posed by Hassolman, Munk & MacDonald [l] as a statistical tool for investi- 
gating nonlinear interacting ocean waves, and Godfrey [2] used th ^ bispectrum 
for analyzing economic time scries. The statistical properties of bispectral 
estimation and variance were reported by Brillinger & Rosenblatt [3] and in- 
cluded an application to the analysis of sunspots. Large sample theoretical 
investigations of bispectral estimates are described in Brillinger ^ 
Rosenblatt [4] , Huber et al. [S] used bi spectra to investigate EEG 
readings I and Roden ^ Bendliner [6] used cross-bispectra to study profiles of 
oceanographic variables such as density and salinity. Various aspects of 
spectral energy transfer in laboratory and atmospheric turbulence have been 
studied using bispectral and cross-bispectral techniques by Lii, Rosenblatt & 
Van Atta [?], by Helland, Lii Rosenblatt ( 8, 9 ], and by Van Atta [lO, ll]. 

Sato Sasaki [12] used bispectral methods to eliminate additive Gaussian noises 
in the reading of laser holograms. A real-time bispectral technique developed 
by Sato, Sasaki & Nakamura [l 3] was used to analyze noise in mechanical gear 
trains, Kim & Powers [ 14] have shown that a bispectral analysis is useful in 
distinguishing between spontaneously excited modes and coupled modes in 
plasma wave studies, and McComas [l3] has reported on theoretical and experi- 
mental bispectra calculations for studies of internal waves in the ocean. 

The purpose of the present study is to review the bias and variance prop- 
erties of the bispectrnm, and consider the practical estimation problem of the 
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general cross-bispectrum for a multivariate time series. The computational 


1 

i 

I 

simplifications which result when the bispectral analysis involves only one or j 

two separate time series are reviewed. A FORTRAN subroutine which imple- \ 

meats these considerations is discussed and a listing of the algorithm is ! 

available from the authors. j 


2. Theory 

We are to estimate the cross -spectrum of three stationary real -valued 
time series defined by X^(i), X 2 (i)t X ^(i); i = 0, ...,T-1 with zero mean. Detailed 
asymptotic results of the k-th order spectra can be found in ( 3] , 

For computational purposes we define the discrete Fourier transform 

T- 1 

F.(X) = ^ X.(t) exp f-i Xt } , j= 1,2,3 (2.1) 

•> t = 0 •* 


and the bispectral periodogram 


R, 2 X X 1 = — JJ F (X 

^ ^ ^ (2tt) T j = 1 J •> 


where 


I 


X . = 0 (modulo 2 n ) 
J 


2 TTS . 

X. with s. integers, 

J ^ J 


( 2 . 2 ) 


If the mean is not removed before taking the Fourier transform, then set 
F.(0) = 0. Now taking the average or expected value of the periodogram 
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(2. 3) 


En 


1,2, 3^^ r ' 3' ^ ^2^ 


where b(X,, X^) is the cross-bispectrum of the joint cumulant of X. ,i = l, 2, 3, 
I s- i 

and E is the expectation operator. Also the asymptotic variance estimate {4^ is 


piven by 


VariRelij ^ ^(Xj, X^, X ^)) 


Var(hnl3j ^ 3 ^^i» ^2’ ^3^^ 




fj(Xi)f2(X2)f3(Xj+ X 2 ) (2.4) 


where f.(X) is the power spectrvim of X. . The asymptotic covariance between 

the real and imaginary parts is an order of magnitude smaller, as is the co- 

variance between different frequencies. To reduce the variance of an estimate 

we can averape B ^ >(^ii ^ from disjoint time blocks (i.e. partition a 

1, 4 ., 3 1 

long stretch of data into disjoint segments or records of equal length) and 
smooth spatially (over frequencies). The spatial smoothing should not be too 
broad otherwise the bias of the estimate would become large. When the bi- 
spectrum is relativelv smooth, the smoothing window shape is not critical. 

F'or most purposes, the uniform wimh^w is adequate and greatly simplifies 
the comput.ition. In the course of computation, it is imporiant to compute an 
estimate of the standard dev'iation of the bi spectrum. If we have m records 
of data and n points being averaged in the frequency domain, then the esti- 
mated standard deviation (S. D. ) at (X j, X.,) for the real part is 


S. D. (X,, X,,) 

1 


, 1 1/2 

VAR Re 13(X ,, X,) 

1 «-> 


m 


(2.5) 


where 


m 


VAR Re B(> ,,>,). -i- Y [Re B. (X , , X ,) - Re BO, , , X ,)]“ (2.6) 
r w 1 1 * Z 1 ' Z " 

1 I 


4 


with 


n 


B.O,. X,) 
I 1 z 




B. 


3<h 




m 



The standard deviation of the imaginary part, Im B, of our estimates is com- 
puted by the analogous method. These estimated standard deviations give us 
a picture of the relative resolvability of bispectrum at each frequency pair 
(> ,, X,) and can be compared with that of asymptotic variance in E^q, (2.4), 

The asymptotic formula Eq, (2, 4) for the variance can serve a fundamental 
role during the design of an experiment. The power spectra, f,(X), are often avail 
able or are readily estimated from preliminary data, 'Fhe asymptotic vari- 
ance formula and preliminary computations of the bispectrum should then be 
used to determine important experimental parameters such as duration of the 
data, number of records, and size of the spatial averaging required to obtain 
a desired degree of statistical resolution in the final bispectral estimate. Tire 
directly estimated standard deviation given in Eq, (2. 5) should be used to verify 
that the degree of statistical resolution has actually been obtained from the data 
analysis. 

3. Symmetries 

The basic symmetry of the cross-bispectrum of real time series is 
obtained by conjugation 

b{> j, ^ b'Vx j, 




(3. 1) 


wher« b* is the complex conjugate of b. Thus, we need only compute the 
right half-plane of Figure 1. When dealing with a single time series, i. e. , 

X = X, = X, , then by various symmetries we only need to compute region 2. 
All other regions can be computed from region 2 using the relationships in 
Table 1. When X^ = X^ = Y, then since 

^XYY^^ l’ ^ 2’ ^ 2^ ” ^XYY^ ' I ’ ^ 2’ ^ 2^ (3. 2) 

we only need to compute regions 1-3, Regions 4, 5, 6 are obtained from 
regions 3, 2, 1, respectively, using Eq. (3. 2) in the appropriate region. 

Note that all these compxitations should be made away from X. = 0 
to avoid introducing a severe bias in the bispectral estimates in the neigh- 
borhood of X . = 0. If estimates are needed at a point D(> ,X X ,) involv- 
ing some > , = 0 then interpolation via a continuity assumption or other 
methods should be made. 

4. Computational Methods 

The fundamental compvitation required for bispectral estimates is 
given by Eq. (2, 2), the triple complex product of the three individvial Fourier 
transforms. While this computation is straightforward, limitations on com- 
puter time, memory, and statistical variance impose severe limitations on 
implementation of the definition of the bispectrum. If record lengths are 
short, then it is a simple matter to compute the triple complex products at 
each pair of frequencies in the right-half-plane. Applications of bispectral 
analyses usually involve long records, and spatial averaging is often 
required in order reduce the variance v f the l)ispectral estimates. 

b 


X 


Increasing the number of records is demanding on computer time and intro- 
duces potential non? tationarities. The spatial averaging can be implemented 
by Slimming over triple prodTicts in a neighborhood of a frequency pair (X 
Ihe summation must be performed record-by-record because the number of 
frequency pairs is large even for moderate record lengths. Computation of a 
large number of b’ spectral estimates is time consuming even on large main- 
frame computers> and it is therefore desirable to compute the estimates over 
an arbitrary set of subregions in the right-half-plane. This approach permits 
taking full advantage of the symmetry relationships discussed in Section 3 . 

In the proposed algorithm the bispectral estimates are computed by 
averaging over rectangular regions which may be located anywhere in the 
right-half-plane and bounded by the Nyquist frequency constraint imposed by 
the discrete Fourier transform. The rectangular regions may be reduced to 
a single point if the variance properties of the data justify the sharp resolu- 
tion. Coordinates for the averaging rectangles are passed to the subroutine 
and the location of the rectangles is arbitrary. This permits easy computa- 
tion of the bispectrum for selected regions when strong interactions are 
expected for isolated frequency combinations. The bispecirum and vari- 
ance estimates are indexed by the averaging rectangle. The bispectrum 
computation is divided into a number of subregions which depend on the loca- 
tion of the averaging rectangle. This is done to eliminate most of the condi- 
tional tests from the inner loops thereby substantially improving the compu- 
tational efficiency. The algorithm requires that the Fourier transform be 
defined out to twice the Nyquist freqx^ency index in order to eliminate 


conditional tests in the first quadrant. Some discrete Fourier transform 
algorithms already contain this symmetrical part and the user is responsible 
for the existence of the extension. A sample computation to extend the trans- 
form beyond the Nyquist frequency is included in the comments of the hi- 
spectral algorithm. 

In the first quadrant the bispectral estimate at a point is given by 

B(i. j) - F^(i) F^Cj) F3(k) (4. n 

where i + j + k = 0, or k = -{i + j). Note that the range of indices i, j, and 
k are assumed to contain 0 and the scale factor l/(2rr)^T has been surpressed; 
the Fortran implementation incorporates a shift by one to avoid the zero index. 

The Ilermitian property of the Fovirier transform of a real time series implies that 

F^dO = F^i-k) (4. 2 ) 

and thus the bispectral estimate in the first quadrant becomes 

B(i.j) = F^(i) F^d) F^C-k) . (4.3) 

As shown in Section 3 bispi'ctral estimates do not need to be made in the 
second and third quadrants. Tn the fourth quadrant i > 0, j < 0, and k may 
be either positive or negative. The dividing line is i = -j. Thus depending 
on the relative magnitudes of i and j, tiie expression for the bispectral 
estimates becomes: For k < 0, 

B(i.j) - F^(i) F^'C-i) F^C-k) , (4.4) 

and for k 0 , 

P(i. l) = F^(i) F*(-j) F^(k) . (4. S) 
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In the implementatJon of these considerations, there is a strip along i = -j 

where k may be of either sign; this is a result of the finite size of the average 

ing rectangles. In the algorithm, the fourth quadrant is divided into three 

separate computation loops, one with k > 0, one with k < 0, and a third with 

k taking on positive and negative values. The point -wise conditional test 

remains in this third loop bvit not in the first and second loops. We feel 

that the small loss of efficiency of including the conditional test to select 

Eq. (4. 4) or (4. 5) in the third lo.)p is worth the simplification gained over 

eliminating even more of the required conditional tests. Tn a global bi- 

spectrum computation only a small percentage of the loops will actually 

include the point -wise conditional test. 

The algorithm does not perform the final scaling for the number of 

records used in the bispectral computation, the conversion to physical units or 
2 

the 1/(2tt) T factor in Eq.(2, 2). The final scaling depends on the particular imple 
mentation of the discrete Fourier transform vised in the computation. 

5. Application 

The proposed method for computation of bispectra is applied to an 
exponentially distributed process. Let X^(il, i = l,...,N be independent, 
identically distrihvited time series with X(i) exponentially distributed with 

parameter a. Then X(i) has a probability function P(X) =— expf-X/a}, 

th k 

X > 0, and the k cumulant of X(i) is C, = (k-1)! Cl . Consider two 

k 


useful test cases: 


(1) If Xj(i) = X^a) = x^(i) then 


V Y ■> » 5^X. = 0 (modulo 2 tt) ; (5.1) 

^ ^ (2tt)^ ^ ' 


(2) If S + application of the running 


mean 


j = - f 


filter to the exponential process with f = (M-l)/2 integer and X^=Xj 


(5.2) 


(X,, X,) = 


C3 Sin(X2M/2) 




(5. 3) 


for -TT < X . s rr. The power spectrum for signal X and X., are 
^ 1 Z 


'x/’' = * • 


and 


(5.4) 


f 1 rSin(XM/2)V 

X,'^' ■ 2rr Lm Sin(X/2)J 


(5.5) 


respectively. 

From Eqs. (5,4) and (5.5) for the power spectra of and X^ the 
theoretical variance of the estimates are calculated from Eq. (2.4) to give 


VAR(b 


^ c,^ 




(5.61 


for case 1 and 


VAR(b ) 


T_ 2 rSin(XM/2) 

4 2 Lm Sin(X/2)J 


.2 


32 TT 


(5.7) 


for case 2. 
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While both cases have been tested, we choose to report on the second 
case because the cross-bispectrum exercises all the loops of the algorithm. 

Case 1 is useful for checking the simple bispectrum computation, but does not 
involve all the possible loops. Case i. was implemented in a test program and 
the cross -bispectrvim was computed using the method described in Section 4. The 
parameters chosen are the transform length T = 256, the number of points included 
in the averaging rectangle n~ 11x11 = 121, and the number of records m = 100. 

In addition, M in Eq. (5, 2) is chosen to be 5, and the exponential process parame- 
ter a= 1. The time series w'as generated and discrete Fourier transformed using 
standard IMSL subroutines. The numerical estimates for the real and imaginary 
parts of the cross-bispectrum together with the standard deviation of the estimates 
are given in Tables 2 and 3, respectively. An extra column is included in Tables 2 
and 3 giving the theoretical values for the real and imaginary parts of the cross- 
bispectrum and the theoretical standard deviations. A single column is sufficient 
since the theoretical values given in Eqs. (5. 3) and (5.7) are constant for 
X 2 = constant. The cross-bispectrum was computed in regions 1-3 of Fig. 1(a) 
and Fig. 1(c). These results are shown in Fig. 2. 

Table 2 shows that the real part of the estimated cross-bispectrum is 
an order of magnitude larger than the estimated standard deviation, and the 
result is considered to be well resolved. Table 3 shows that the estimated 
imaginary part is effectively zero since it is comparable to or smaller than 
the estimated standard deviation. A comparison of estimated and theoretical 
values for the real part in Table 2 shows that the differences are small relative 
to the staiulard deviation whereas the estimated and theoretical standard 

i 



deviations are comparable. The comparison between numerical and theoretical 
cross-bisoectra has been statistics', in nature and clearly exact agreement 
should not be expected. The accuracy of the algorithm was tested by a deter- 
ministic sine-cosi..e signal, as wall as the application just discussed. The 
comparison obtained between the numerical and theoretical cross-bispectra 
provide strong support for the accuracy of the algorithm in estimating both 
the mean and standard deviation of bispectra and cross-bispectra. 

6, Concluding Remarks 

We have described a method for computing the bispectral and cross- 
bispectral estimates of up to three separate time series at any point in the 
right-half-plane. Proper application of the algorithm involves: 

(i) The computational domain must be chosen to minimize the computation 
time required by observing the symmetry relationships. 

(ii) The averaging rectangles should be of optimum size to both keep the 
ratio of mean to standard deviation as large as possible without making the 
rectangles so large as to lose the spatial (frequency) resolution of important 
features in a particular application. If there is no theoretical guidance avail- 
able in a particular application, then the user must proceed in a trial-and 
error technique to find the optimum computational parameters. 

A source list.ng of the computer subroutine which implements the pro- 
posed method is available from the authors. 


12 


Acknowledgment j 

The authors wish to thaiJ< Professor M, Rosenblatt and Mr. R. A. 
Stanford for helpful discussions. This work was supported by NASA Grant 
NSG 237 6. KNH was partially supported by National Science Foundation Grant 
ENG78 -25080. 


13 


References 


1. Hasselmann, K. , Munk, W. , and MacDonald, G. Bispectra of Ocean 
Waves. In Time Series Analysis (ed. M, Rosenblatt) (196 3), Wiley, 
125-139. 

2. Godfrey, M. D. An exploratory study of the bispectrum of economic 
time series. Appl, Statist. (London) 1 4(1965), 48-69. 

3. Brillinger, D. R, , and Rosenblatt, M. Computation and Interpretation 
of kth order spectra. In Spectral Analysis of Time Series (ed. B. Harris) 
(1967a), Wiley, 189-2 32. 

4. Brillinger, D. R. , and Rosenblatt, M. Asymptotic theory of estimates 

of kth order spectra. In Spectral Analysis of Time Series (ed. B. Harris) 
(1967b), Wiley, 

5. Huber, P. J., Kleiner, B., Gasser, Th. , and Dumermuth, G. Statistical 
methods for investigating phase relations in stationary stochastic proc- 
esses. IEEE Trans, Audio Electroacoust. 19 (1971), 78-86. 

6. Roden, G. I. , and Bendiner, D. J. Bispectra and cross-bispectra of 
temperature, salinity, sound velocity and density fluctuations "th depth 
off Northeastern Japan, J. Phys. Ocean 3, 3 (1973), 308-317. 

7. Lii, K. S. , Rosenblatt, M. , and Van Atta, C. W. Bispectral measure- 
ments in turbulence. J. Fluid Mech . 1 (1976), 45-62, 

8. Helland, K. N. , Lii, K. S. , and Rosenblatt, M. Bispectra of atmos- 
spheric and wind tunnel turbulence. In Applications of Statistics 
(ed. P. R. Krishnaiah) (1977), North- Holland Pub. , 223-248. 


14 



9. Holland, K. N, , Lii, K, S. , and Rosenblatt, M. Bispectra and energy 
transfer in grid- generated turbulence. In Development in Statistics 
(ed. P. R, Krishnaiah) 2 (1979), 123-155. 

10. Van Atta, C, W. Bispcctral measurements in turbulence computations. 
Proc. 6th International Conference on Numerical Methods in Fluid 
Dynamics, In Lecture Notes in Physics (1978), Springer- Verlag, 
530-536. 

11. Van Atta, C. W. Inertial range bispectra in turbulence. Phys. Fluids 

8 (1979), 1440-1442. 

1 2. Sato, T. , and Sasaki, K. Bispectral holography. J, Acoust. Soc . 

Am. 2 (1977), 404-408. 

13. Sato, T. , Sasaki, K. , and Nakamura, Y. Real-time bispectral analysis 
of gear noise and its application to contactless diagnosis. J. Acoust . 
Soc. Am . 2 (1977), 382-387. 

14. Kim, Y. C., and Powers, E. J. Digital bispectral analysis of self- 
excited fluctuation spectra. Phys. Fluids 21 , 8 (1978), 1452-1453. 

15. McComas, C, H. Bispectra of internal waves. Woods Hole Technical 
Note 02543 (1078), Woods Hole, Mass. 


15 


Table 1. Bispectral symmetries 


b(X^, X^) in 


Region 1 

3 

4 

5 

6 


Corresponding bispectral coordinates 
in Region 2 


b(X 2- Xj) 


b"^(Xi+ X^. 


b (.X2.X^ + X2) 


bCX^.-X^-X^) 


M-X^-X^. X^) 


U. 


<\J 

r- <M 

<Vi 

f-i M 

M 

(NJ CO 

^ rO 

(M CO 

o o 

o o 

O O 

o o 

O O 
1 1 

o o 
1 1 

o o 
1 1 

O O 
1 1 

U U 

< 1 
u u 

U U 

U U 

U U 

U U 

u u 

U M 

tn vO 

m M 

CM t^ 


O 00 

O 00 

fO in 

f*« fO 

fM fO 

Ov Tf 

o ^ 

cn 

nj O' 

o 

ro CO 

in 00 

vD ro 

00 

m 

in 

(VJ 

ir> in 

O nj 

•-H ^ 

ro 

<M (NJ 

Tf in 

CO fO 

t— • r- « 

00 r- 

.-t (T' 

no no 


nj 

ro 


no 

t-^ 

no 

no 

CO 


no 

•-I 

ro 

o 

O 

o 

o 

o 

O 

o 

o 

o 

o 

o 

1 

o 

1 

U 

u 

u 

u 

W 

w 

u 

u 

u 

u 

u 

u 

no 

CO 


no 

vT* 

in 

O 

O 

no 

in 

no 


'•T 


CO 

o 

vC 

o 

r— 


00 


o 

O 

in 


o 

CO 

CO 



ro 

a- 

ro 

cn 

nC 


nj 

•—I 

• 

r- • 

« 

— 


<r—* 

ro 

no 

no 

• 

• 

• 


no 

ro 


no 

no 

CO 

r— • 

no 


no 


ro 

r-i 

no 

no 

CO 

o 

O 

o 

O 

o 

o 

o 

o 

o 

o 

O 

o 

o 

1 

o 

1 

o 

1 

o 

1 

1 

U\ 

U 

u 

u 

w 

u 

u 

u 

u 

w 


w 

u 

u 

w 

u 

cn 

r'~ 


O 

O' 

rO 

ro 

in 

00 

r- 


O 

o 

CO 

O 

in 


ro 


cr* 

fO 

no 


in 

r-* 

o 

cr^ 

00 

o 

O 


o 

vO 

»—* 

ro 

ro 


in 

CO' 

no 

in 

ro 

O' 

cn 

ro 

O 


n* 

cn 

00 

• 

O' 

*—• 

no 

no 

no 


ro 

rO 

ro 

r-< 

* 

f-^ 

f' 

vO 



CO 

*-• 

no 

no 

CO 

— CO 

•— < 

no 

*— * 

ro 

»— * 

ro 

o 

o 

o 

o 

o 

O 

o o 

o 

o 

o 

o 

o 

1 

o 

u 

w 

u 

u 

u 

M 

u u 

u 


u 

w 

u 

u 

vC 

r- 



vO 

n- 

ro 0> 

CO 


TT 

00 


no 

CO 

cc 

00 

CO 



in tT 

00 


O' 


00 

O" 

O 

n- 

nj 

ro 

vO 

r— » 

C' 


or 

ro 

r- 

CO 

in 


« 

r- • 

• 

• 


CO 

no no 
• • 

in 


• 

CO 

• 

r-» 




no 

ro 


CO 


rsj 

no 

ro 

f— • 

no 

•-H 

ro 


ro 


o 

O 

c 

O 

o 

o 

o 

O 

o 

o 

o 

o 

1 

o 

1 

o 

1 

O' 

o 

1 

W 

U 

W 

w 

W 

W 

W 


W 

W 

U 


W 

U 

ro 

no 

Vk.' 

hO 

ro 

<X) 

in 

ro 


oa 

ro 

in 

in 


t- 

• 

TT 

nl 


ou 

n- 

r- 

00 

o 


no 

ro 

O' 

n- 

o 

o 


no 

c 

vC 

nj 


ro 

00 

O' 

no 

r-» 

oc 

00 

a> 


f— < 

no 


00 


r— ' 


no 

no 

no 

in 

CO 

ro 

rO 








f— » 



OC 


oc 




no 

O' 

0 

ro 

o 

cc 


ro 

o 

ro 


r- 

• 

• 

d 

d 

d 

d 

• 

o 

d 

d 



04 

03 

03 

04 

02 

CM 

O 

-02 

CO 

o 

1 

CO 

o 

1 

CO 

o 

1 

CO 

o 

1 

1 

u 

OD 

O r-* 
vO 
00 

• 

1 

w 

00 

sD 

in 

• 

0 

.6450E- 

0 

.6657E- 

0 

. 1530E- 

U 

o 

m 

CM 

• 

M 

O ^ 

CJ 

• 

0 

.7012E 

W 

vC 

• 

M 

CJ 

m 

t 

u 

Ift 

• 

04 

03 

03 

03 

CO 
O O 

rO CO 
O O 

-03 

-02 

CO CM 

o o 

CM CM 

o o 
1 1 

-03 

-03 

ro 
O O 
1 1 

-02 

-03 

-03 

-03 

.6328E- 
. 1833E- 

1 1 

u u 
00 

in 

in 

v£) vO 

• « 

1 

U U 

00 o 
vO in 

CM 
m 00 

• * 

u u 

so r- 
O nO 

•- 

CO 

ft ft 

.607 2E- 
. 1666 E- 

-.787 3E. 
. 2328E- 

M U 

CM 

CM 
Cv o 
r-i CM 

• • 

w w 

CO CM 
CM M* 
vO - 
sD C ^ 

« • 

1 

w w 

00 O' 
o in 

CO vO 

• » 

w w 

Tf CM 
sC nT 
CO 
r-l sO 

« • 

1 

U K 

QC 

r- c 
in CM 

fVl r- 

ft 

03 

03 

CM 
O O 

m CO 

o o 

.03 
• 03 

-02 

-02 

CO CM 

o o 

CM CM 
O O 

1 1 

CM CO 
O O 
1 1 

-03 

-03 

-03 

-03 


2077E- 

1928E- 

5656E- 

lllOE- 

1389E. 

7612E- 

2345E- 

2377E- 

w w 

r- 
in CM 

o r- 

r-i »** 

U U 

M* r- 
o 

u u 

vO C30 
CM 

CM 
CM CM 

u u 

nO CO 
so 

in in 

9^ 00 

3593E 

5141E 

367 IE 
266 5E 

ft 

• • 

* • 
1 

ft ft 

t 

1 

ft ft 

ft ft 


1 

1 

1 

o c 

CO 
O O 

03 

03 

CO CO 

o o 

CO CO 

O o 

CM 

o o 

CO CM 
O O 

-02 

-02 

tT CO 

o o 

t 1 

-03 

-03 



8761E- 

1949E- 

1 1 

w w 

00 
00 
1— 1 m 

w w 

o CM 
O O' 

r-« 

m vD 

w w 

cr CM 
m vo 
CM f>- 
rO CM 

3564E- 

1650E- 

8030E 
247 7 E 

W W 

M3 CO 

O o 

in CM 
CM CM 

.7 447E 
.8888E 

2005E 
2287 E 


ft 

# • 

ft • 

1 

ft ft 

1 

ft ft 

1 

ft ft 

1 

1 

1 

1 


o o 

o c 

03 

03 

03 

03 

CM CO 

o o 

Tf CO 

o o 

CO CM 
O O 
1 1 

CM 

o c 

n* CM 

0*3 o 

t 1 

-03 

-03 




i 1 

w w 

O 00 
iTi O 

lO »— • 

,-4 rsj 

3264E- 

6321E- 

w w 

CM 00 
CO O 

Tf o 

i— « 00 

W W 

r-i 

^ cn 

m f>- 
CM 

w w 

o •-* 
in CO 
O' r- 

CO »-i 

W W 

CM nO 
00 

CM t\i 

5318E 

3122E 

M Ui 
in 
cr* vO 
sO c 

CM ro 


« % 

ft 

* « 

ft • 

ft ft 

1 

ft ft 

1 

• • 




o o 

0 

0 

o c 

04 

03 

03 

03 

• 03 
03 

-03 

-03 

-02 

-02 

-02 

-02 

-03 

-03 





6454E- 

1667E- 

i 1 

u u 

O' ^ 

7749E- 

5708E- 

3441E- 

2351E- 

3163E- 
. 1578E 

w w 

00 CM 
O 00 
O CM 
CM 

5043E 

8090E 

ft ft 

ft ft 

• • 

ft 




Figure Captions 


b^spectral domains for discretely sampled data: 

(a) Xj + + X j •- Oj (b) X j + X 2 + ^ j = -2 tt; 

(c) ^ + ^2 ^ ^ 3 ~ 

Comparison r ' estimated and theoretical cross-bispectra, 
(a) Theoretical real part; (b) estimated real part; 

(ci estimated imaginary part. 





